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1. Motivation and Overview 

The motivation to study excited-state contamination in nucleon correlation functions using 
chiral perturbation theory comes directly from lattice QCD. As a result, our primary concern is 
with the Euclidean time dependence of correlation functions, such as two-point and three-point 
functions that arc relevant to extract the nucleon mass and current matrix elements using standard 
lattice QCD techniques. For a rundown of such techniques relevant to baryon physics, see [1], In 
the case of the nucleon two-point function, large Euclidean time separation is used as a filter for 
the ground-state nucleon. From the zero three-momentum projected correlator, G(x), one typically 
looks at the effective mass 

M eff (T) = -log + T =“ M n + I Z'\ 2 e- AEx + • • •, (1.1) 

G{ T) 

where a is the lattice spacing, Mn is the nucleon mass, and A E is the energy splitting with the near¬ 
est excited state, A E = E' — Mn, with Z' | 2 as the ratio of the excited-state to ground-state coupling 
of the chosen lattice nucleon interpolating operator. A flat effective mass indicates that the ground 
state saturates the correlation function, and allows for a reliable extraction of the nucleon mass. To 
reduce uncertainty surrounding saturation of the ground state, one should go to suitably long Eu¬ 
clidean time separations, however, lattice QCD correlation functions are determined stochastically. 
For baryon correlation functions, it is well known [2] that the signal in Eq. (1.1) degrades expo¬ 
nentially compared with the statistical noise. Without an increase in statistics, one can attempt to 
minimize |Z'| 2 through the construction of better nucleon interpolating operators. Such operators 
couple less strongly to excited states; and, to this end, we investigate the nature of pion-nucleon 
excited states using chiral perturbation theory. One should note that the signal-to-noise problem 
is more restrictive in the case of three-point functions, due to the necessity of two Euclidean time 
separations that should be ideally long. 

Here we address the pion-nucleon (and, to a lesser extent, pion-delta) contributions to nucleon 
correlation functions. The un-amputated correlators are computed as a function of Euclidean time 
using chiral perturbation theory. This is a summary of previous work, namely: the computation of 
excited-state contamination using a continuum of excited states was pursued in [3], while this com¬ 
putation was recently revisited for the case of discrete pion-nucleon and pion-delta states available 
on a periodic lattice [4], We detail the computation of the nucleon two-point correlation function, 
because it contains universal couplings to pion-nucleon states. We also investigate the excited-state 
contamination present in the nucleon three-point function of the isovector axial-vector current. The 
forward matrix element of this current is parameterized by the nucleon axial charge, gA■ We find the 
sign of excited-state contamination in gA is inconsistent with underestimation of the axial charge. 
An outlook to future work is given in the summary. 

2. Nucleon Operators 

To begin, we note that within chiral perturbation theory the nucleon is described by a local 
operator, N(x). Of course locality is only defined relative to some scale, and this scale is set 
by the inverse of the chiral symmetry breaking scale A x = 2\/2nf, where / is the pion decay 
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constant and / = 132 Me V. The long-distance nature of the pion means that observables in chiral 
perturbation theory are computed as long-range pion contributions plus short-range physics, see 
Fig. 1. In the effective field theory, we know that the nucleon operator transforms as a doublet 
under the unbroken isospin symmetry, SU (2)y. It is, however, unknown to which chiral multiplet(s) 
the nucleon belongs in the chiral limit. In fact, the transformation property of the nucleon field 
under SU(2)l x SU(2)r is not an ingredient of chiral perturbation theory. It is both customary 
and convenient to have the nucleon field transform as N(x) —> U (l.R. q(x)'jN(x) under L. R £ 

SU(2)l x SU{2)r, where £ is the square-root of the coset field, q = \/E, and transforms as ^ M 
Ut,R' = LqU'. This freedom stems from the ambiguity in resolving a nucleon and any number of 
pion states in the chiral limit. Using this simple transformation, the pion-nucleon (and pion-delta) 
interactions are then easily constrained by the pattern of spontaneous and explicit chiral symmetry 
breaking of QCD. 



Figure 1: Graphical depiction of the nucleon operator in chiral perturbation theory, as well as a lattice QCD 
interpolating operator for the nucleon. In the former, the size of the short-distance core is set by Ar , while 
in the latter, the operator smearing radius R@ sets the size. 


In lattice QCD, the situation is actually quite analogous. One studies properties of the nucleon 
using a lattice interpolating operator, &n{x). To overlap with the nucleon, this operator must have 
the correct isospin quantum numbers. Unlike chiral perturbation theory, however, lattice interpo¬ 
lating operators belong to a given chiral multiplet. The chiral properties of baryon operators have 
been studied completely [7]. We will use the most common of such operators 

& N ~ q (q T C'Y5% 2 q ), ( 2 . 1 ) 


where t are the isospin matrices, and q is an isodoublet of quark fields, q 



. This operator 


transforms as (2l, 1#) © (1^,2^) under chiral transformations. Typically nucleon operators are 
chosen to maximize overlap with the ground-state nucleon, which is characterized by the overlap 
factor 


Ztf = 


( 2 . 2 ) 


There are also excited-state couplings of the operator, for example Z k n = (0| (nN) p ), where the 

p denotes a relative p- wave. One way to maximize overlap with the ground state is to smear the 
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operator over some size Rtf. It turns out that the relative contribution of the pion-nucleon excited 
states compared to the ground state, \Z n ^\ 2 /\Z^\ 2 , is constrained by chiral dynamics provided the 
operator smearing radius is small compared to the pion Compton wavelength, Rtf <C m 1 [5]. A 
similar idea was put forth for gradient flow observables by considering the flow time as a scale in 
chiral perturbation theory, see [6]. 

The above discussion implies the procedure to compute the chiral contamination in nucleon 
correlation functions. For a given lattice interpolating operator, &n{x), we assume the hierarchy of 
scales R@ V A x 1 < m n 1 , so that the long-range physics remains dominated by pions. We 
then employ chiral perturbation theory for the lattice operator, by mapping the lattice interpolating 
field into effective field theory operators, Ntf = Ntf(N. q). This mapping is also required in the case 
of chiral perturbation theory for moments of the nucleon distribution amplitudes [8]. In that case, 
of course, the nucleon operators are strictly local. 1 We then use the effective field theory nucleon 
operators to compute correlation functions. As an example, nucleon chiral doublet operators take 
the form 

N {2l ,U) ~ ^iV + tf(V), and N{W~\Vn + 0{M^). (2.3) 

This leads to the follow heavy nucleon effective field theory operator for the lattice interpolator 

N<2 L ,1 R M1 L , 2 R ) = \ztf(^ + ? )N = Ztf(l-^ + ---^N. (2.4) 

From this expression, we can see that the relative coupling of pions to the nucleon is fixed for 
the given lattice operator because of that operator’s chiral transformation properties. For a given 
nucleon correlation function, there will be pion loop diagrams that contribute in addition to the 
pions generated from the chiral structure of the lattice operator. 


3. Chiral Contamination in the Nucleon Two-Point Function 


Now we turn our attention to the determination of the nucleon two-point function using chiral 
perturbation theory. This two-point function has a matrix element definition 

G( T) = I(0| 0 N (x, t) 4(0,0)|0), (3.1) 

X 

where the sum over all lattice sites x projects the correlation function onto vanishing three mo¬ 
mentum. Quite generally the two-point function has a spectral representation with positive spectral 
weights representing the probability to find excited states in the correlator. For our purposes, we 
write the spectral representation in the form 


G{ t) = \Ztf\ 2 e~ MNT 


i + Ei 

n=f^ 0 


-AE„t 


(3.2) 


where the excited-state couplings |Z„| 2 are defined relative to the ground-state nucleon, as is the 
energy A E n , which has the form A E„ = E n — Eq, where Eq = M # is the ground-state energy, which 

'Chiral perturbation theory calculations of excited-state contamination have been performed for two-point correla¬ 
tion functions of local axial-vector and pseudoscalar bilinear operators in vacuum [9]. 
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is just the nucleon mass. Chiral perturbation theory makes a prediction for the relative contribution, 
|Z„| 2 , from the various low-lying excited states. Our primary concern is with intermediate pion- 
nucleon states, see Fig. 2. 




Figure 2: Feynman diagrams depicting excited-state contamination in the nucleon two-point function. A 
single pion shown on the left, while an intermediate state with multiple pions is shown on the right. 


The spectral density has previously been determined in chiral perturbation theory [3]. In that 
work, a continuum of pion-nucleon states was employed, for which one has the relation £„ |Z„| 2 —> 
/“ dE p n u{E). The spectral density has the form 


PME) 


6 gj for | 3 

(47T/) 2 fc2+ m 2 


(3-3) 


which has a number of features. First is the non-analyticity at threshold, such non-analyticity can 
only arise from a continuum of states. Second is the behavior near threshold, |k^| 3 , which arises 
because of the two-body phase space ~ \k K \, and the requirement that the nucleon and pion be in 
a relative p- wave, ~ k 2 . By contrast, on a periodic spatial lattice of length L in each direction, the 
allowed pion momentum modes are restricted to k n = =j-n. For such modes, we have the weights [4] 


-itN\ 


3 sl kj 

4 f 2Ll [kl + ml p/ 2 ’ 


(3-4) 


which are analytic at threshold; because, there, we no longer have an accumulation of states. Thus 
at threshold, the behavior of the two-body relative weights is determined solely by the requirements 
of angular momentum. It should be noted that the two-point function of the nucleon has been 
determined using relativistic baryon chiral perturbation theory [5]. Expanding these results in 
m n /M N 0.15, they reduce to those given in Eq. 3.4. 

The effective mass of the two-point correlation function in Eq. 3.2 is plotted in Fig. 3. For the 
relative spectral weights, we use those derived from chiral pertubation theory, Eq. 3.4. To make 
the plot, we use the physical value of the pion mass, and choose a lattice spacing of a = 0.1 f m on 
a lattice of size L = 48 a. From the plot, we see that pion-nucleon excited states do not contribute 
substantially to the two-point correlation function. In fact, the suppression of pion-nucleon states, 
and more generally (pion)"-nucleon states (see Fig. 2) can be argued without any reference to chiral 
perturbation theory. The reasoning stems from simple phase-space considerations near threshold. 
For a continuum of states, the spectral density near threshold is determined entirely by the available 
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Figure 3: Effective mass of the nucleon two-point function computed with chiral perturbation theory. Pion- 
nucleon excited-state contamination is shown at the few-percent level for typical values of lattice parameters. 
For times such that x/a>5, the correlator is saturated by pion-nucleon states with momentum mode numbers 
satisfying ViP < 4. 
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phase space, and the orbital angular momentum t of the multi-particle state. As a result, we have 2 

p [K y, N (E) « h 3n+2e ~ 2 . (3.5) 

Indeed there is an analogous result for the spectral weights of such n -pion states available on a spa¬ 
tial torus. The suppression of such multi-particle states is comparatively less than for a continuum 
of states, as the weights take the form 3 

\Z (k) » n \ 2 ~ VE^E^ n+2i ~\ (3.6) 

Chiral perturbation theory implies that these multi-pion states are further suppressed by factors of 

C si/fL) 2n ■ 

Phase space arguments near threshold require non-relativistic pion states, i.e. we must satisfy 

p 

E — Ej h = H-< 1, which translates into the condition m K L 2> 1. Suppression of the weights 

near threshold, however, can compete with exponential enhancement of excited states in the corre¬ 
lation function. This is because enlarging the volume to approach threshold, also reduces the energy 
of the multi-particle state. As the gap is lowered, the exponential factor e Af - T increases. This may 
or may not be an issue depending on how one approaches the chiral and infinite volume limits. In 
the case of fixed pion mass and increasing lattice volume, there is no issue because m n L increases, 
and the decrease in the gap will be negligible compared to m n . On the other hand, lowering the 
pion mass at fixed lattice size is problematic. In this case, both m K L and the gap arc decreasing. 
Both conditions lead to enhancement of multi-particle states. In this situation, however, one should 
generally be worried about finite volume effects to begin with. 

2 Here we define n to be the number of pions, whereas in [4], we use the integer N to denote the number of particles 
in the intermediate state, which is thus n + 1 for the present case. 

3 One interesting feature is that there is no threshold suppression from phase space of JTiV-states in a relative s-wave 
on a torus. 
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4. Axial Current in the Nucleon 

Given the above computation of the two-point function, we now turn our attention to analogous 
matters in the case of nucleon matrix elements of the iso-vector axial-vector current. For simplicity, 
we consider the forward matrix element which is directly proportional to the nucleon axial charge, 
gA- The un-amputated, three-point function of the axial current between nucleon states is defined 
by 

(0|^(t)/+(04(0)|0> = {w^hW)^aM, (4.1) 

where the iso-vector axial-vector current is given by /+ = qju 75 t 1 q , and we have suppressed 
the projection onto vanishing three-momentum. The rest-frame spinors are denoted by \ff, and ,S U 
is the covariant spin operator. The axial current matrix element depends on two Euclidean time 
separations: the source-to-sink separation z, and the source-to-current insertion separation, t. 



Figure 4: Feynman diagrams depicting the dominant contribution to excited-state contamination in the 
nucleon three-point function of the axial current. 


The standard lattice QCD method to compute the axial charge relies on taking the ratio of 
three-point to two-point correlation functions. In chiral perturbation theory, we define the analo¬ 
gous ratio 

G A {z,t) = y A {z,t)/G{z). (4.2) 

The time-dependence of this ratio is generally given by an expression of the form 

G A (z,t)=gA+ £ fl (4.3) 

n,m^( 0 , 0 ) 

where g„ m arc the axial transition couplings between the n -th and m-th excited states. And the 
excited state energies, AE n are measured relative to the nucleon mass. For Euclidean time separa¬ 
tions t and z — t that are not ideally long, the leading excited-state contamination in the three-point 
function arises from N-to-nN transition couplings, as shown in Fig. 4. 

In Fig. 5, we plot the chiral perturbation theory evaluation of the three-point function ratio 
Ga("M). Notice that in the dual limit of z l$> t and t > 1, we have G A (z,t) -> g A - What is of 

interest concerns how the three-point function ratio approaches the axial charge, i.e. from above 

or below. Unlike the two-point function, the couplings g n>m are real valued, and therefore not 
necessarily positive. Looking at Fig. 5, we see that pion-nucleon contamination drives the axial 
current correlation function upwards. In the plot, we show the axial current correlator ratio as a 
function of the current insertion time t for three values of the source-sink separation, z/a = 10 , 
12, and 14, with a = 0.1 fm taken as a typical value for the lattice spacing, on a lattice of size 
L = 48n. Because the pion-nucleon excited states drive the axial current correlator upwards, these 


7 






Excited-state contamination £ nucleon correlators 


Brian C. Tiburzi 


1.10 

. 

. 

< 1.08 
W) 

■ © A 

A © ■ 

^ 1.06 
■4—> 

■ © A 

A © ■ 

F- 1.04 



< 

O 1.02 

1.00 

. © A 

" ■ f| 

_i_1_i_i_i_1_i_i_i_1_i_i_i_ 

■©► 

■ © ► 

■ © 

■ 


-6 -4 -2 0 2 4 6 


▲ r/a = 10 
© r/a = 12 
■ r/a = 14 



T 

2 


)/a 


Figure 5: Excited-state contamination in the axial-vector three-point function of the nucleon. The con¬ 
tamination function Ga (x.t) is plotted as a function of the current insertion time t (measured relative to 
the midpoint between source and sink), for three values of the source-sink separation, x/a = 10, 12, and 
14. Excited states with pions tend to drive the axial correlator upwards, which could potentially lead to an 
overestimation of gA of a few percent at most. 


contributions could lead to an overestimation of the axial charge if not properly accounted for. This 
trend is opposite lattice QCD data for the axial charge, which overwhelmingly approach a plateau 
from below. For a summary of lattice QCD results for the nucleon, see [10, 11]. 

5. Summary and Outlook 

Above we see how chiral perturbation theory can be used to address pion-nucleon and pion- 
delta contributions to nucleon correlation functions. The applicability of chiral perturbation theory 
for the lattice nucleon interpolating fields requires that the smearing radius of the lattice operator 
is considerably smaller than the Compton wavelength of the pion, R& <C m K 1 . In this regime, the 
lattice interpolating field can be systematically mapped into local operators in chiral perturbation 
theory. This mapping gives one the ability to compute excited-state contamination to lattice nucleon 
correlation functions. In particular, the perturbative pion-nucleon and pion-delta contamination. 

In considering the nucleon two-point function in a general context, we find that phase-space 
arguments strongly disfavor large pion-nucleon contamination. In chiral perturbation theory, there 
is additional suppression due to the strength of pion-nucleon interactions, ~ g A /(fL) 2 . For weak 
pion-nucleon interactions, these results suggest that excited-state contributions to nucleon two- 
point functions would correspond to excited-state hadrons (in the case where they are bound states). 
Consideration of the axial-vector current matrix elements in the nucleon, we see that the dominant 
excited-state contamination arises from nucleon to pion-nucleon intermediate states. For insuffi¬ 
cient time separation between source, sink, and current insertion time, such excited-state contam¬ 
ination drives the axial correlation functions upward. This would lead to an overestimation of the 
nucleon axial charge, g&, if such contributions are not properly removed. On the basis of the per¬ 
turbative pion-nucleon interactions in chiral perturbation theory, we conclude that underestimation 
of the nucleon axial charge is not likely due to pion-nucleon excited-state contamination. 
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Looking to future work, it would be interesting to compute the excited-state contributions to 
other observables determined from lattice three-point functions, such as the nucleon scalar and 
tensor charges, quark momentum fraction in the nucleon, etc. These results, at minimum, would 
provide an indication of the sign of excited-state contamination in these quantities, and this can 
be compared with the actual behavior of the lattice correlation functions. Furthermore, lattice 
QCD studies of correlation functions using well-localized nucleon interpolating fields would allow 
confrontation with these predictions from chiral perturbation theory. For the case of resonances, 
our work shows the necessity of including pion-nucleon operators. This is, by now, not at all 
surprising to lattice practitioners. The case of resonances is additionally interesting due to the 
non-perturbative nature of the pion-nucleon interactions. In this respect, it would be beneficial 
to understand nucleon, pion-nucleon, and resonance contributions to correlation functions in the 
context of unitarized chiral perturbation theory models. 
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